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SUMMARY 

Few  nonintrusive  techniques  are  available  for  particle  measurements  in  the  submicron 
size  range  (<  0.1  ym  diameter),  yet  measurement  of  these  particles  is  basic  to  an 
understanding  of  important  processes  in  combustion,  such  as  soot  formation  and  oxi¬ 
dation.  The  objective  of  the  present  research  is  the  development  and  application  of 
a  technique  for  measurement  of  individual  submicron  particles  in  a  gas  stream.  The 
approach  is  to  measure  the  inertial  relaxation  time  of  individual  particles  in 
Brownian  motion,  by  statistical  analysis  of  the  time-resolved  (100  MHz)  heterodyne 
signal  obtained  in  an  interferometer  system  resembling  a  more  conventional  laser 
|  Doppler  velocimeter.  Progress  to  date  has  included  the  development  of  an  optical 
I  system  and  fast  data  acquisition  system  for  conducting  the  experimental  studies.  A 
detailed  Monte  Carlo  numerical  simulation  of  the  Brownian  motion  of  isolated  particles 
j  has  been  developed  and  applied  to  the  calculation  of  the  time  dependent  signals 
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19.  Abstract  (Continued) 

V  ' 

expected  in  the  experiment.  These  simulated  signals  have  been  analyzed  for  correla¬ 
tions  between  the  particle  size  and  statistical  properties  of  the  signal.  The  mean 
time  between  signal  extrema  has  been  found  to  be  proportional  to  the  square  root  of  the 
product  of  the  particle  relaxation  time,  a  quantity  which  is  proportional  to  the 
particle  diameter,  and  the  data  sampling  interval. 

Continuing  research  on  this  program  will  be  initially  focused  on  validation  of 
the  experimental  approach,  by  statistical  analysis  of  signals  obtained  with  the 
interferometer  for  particles  of  known  size.  Following  the  demonstration  of  this 
approach,  the  technique  will  be  applied  to  the  study  of  the  size  evolution  of 
combustion  particulates. 
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1.  RESEARCH  OBJECTIVES  (STATEMENT  OF  WORK) 

The  research  objectives  for  this  program  are  outlined  In  the  Statement  of 
Work,  which  Includes  the  following  tasks: 

Task  I  -  The  contractor  shall  construct  a  laser  diagnostic  apparatus  to  be 
used  for  submicron  particle  size  analysis  by  measurement  of 
Brownian  motion.  The  development  of  this  apparatus  shall  Include: 

a)  Construction  and  assembly  of  relevant  optical  components; 

b)  Adaptation  of  suitable  data  acquisition  and  recording 
apparatus;  and 

c)  Development  of  appropriate  software  for  analysis  of  signals 
obtained  from  the  Brownian  motion  Instrument. 

Task  II  -  Using  the  optical  apparatus  and  data  recording  and  analysis  tools 
developed  In  Task  I,  the  contractor  shall  investigate  the 
development  of  optical  heterodyne  techniques  to  be  used  in 
measurement  of  the  Brownian  motion  of  Individual  particles. 
Specifically,  techniques  are  required  which  permit  measurements  of 
motion  where  the  particle  velocity  changes  in  distances  shorter 
than  one  light  wavelength. 

Task  III  -  Using  the  Instrumentation  and  techniques  developed  in  Tasks  I  and 
II,  the  contractor  shall  analyze  the  statistics  of  particle  motion 
in  the  gas  to  investigate  the  particle  mass,  friction  coefficient, 
and  equilibrium  temperature.  These  studies  shall  be  conducted 
under  carefully  controlled  flow  conditions,  initially  using 
standard  size  particles.  Subsequent  development  of  the  technique 
may  utilize  a  well-characterized  flat-flame  burner. 


Task  IV  -  The  contractor  shall  investigate  the  measurement  of  the  refractive 
index  of  particles  in  the  Rayleigh  scattering  regime  by  analysis  of 
the  scattered  light  intensity  and  particle  mass  (size)  obtained 
using  the  Brownian  motion  sensor. 

Task  V  -  By  application  of  the  Brownian  motion  sensor  in  a  flat-flame 
burner,  the  contractor  shall  assess  the  application  of  this 
technique  for  in-sltu  sizing  of  submicron  particles  in  combustion 
s  treams . 
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2.  RESEARCH  STATUS 


The  research  effort  during  the  first  year  of  the  program  has  focused  on 
development  of  the  theory  of  single  particle  sizing  by  measurement  of  Brownian 
motion.  The  approach  has  been  to  use  a  numerical  simulation  of  the  Brownian 
motion  of  individual  particles  to  predict  the  behavior  (temporal 
characteristics)  of  the  signals  which  would  be  measured  by  an  experimental 
device.  This  approach  has  proven  to  be  extremely  useful,  both  in  devising 
methods  of  analysis  of  the  signal  statistics  for  determination  of  particle 
size  and  in  guiding  the  design  of  an  experimental  apparatus  (including  data 
acquisition  system)  for  testing  the  measurement  concept.  This  model  of  the 
experimental  technique  is  described  in  detail  in  the  Appendix  to  this  report, 
and  the  important  results  from  the  simulation  as  they  bear  on  the  experiments 
are  summarized  below. 

An  experimental  optical  system  and  a  data  acquisition  system  for  the 
Brownian  motion  measurements  were  assembled  during  the  first  year  of  the 
program.  These  systems  are  described  below.  Finally,  the  research  plan  for 
the  second  year  of  the  program  is  outlined  in  this  section. 

2.1  Theory  of  Measurement  Technique 

Much  of  the  first-year  effort  was  concentrated  on  development  of  the 
theory  of  the  Brownian  motion  measurement  technique.  This  effort  included  the 
development  of  a  detailed  Monte  Carlo  simulation  of  the  experimental  system. 
The  important  feature  of  this  model  is  that  it  simulates  time  dependent 
signals  arising  from  particles  of  different  size  and  facilitates  the 
development  of  approaches  to  signal  analysis  in  order  to  discover  correlations 
between  the  particle  size  and  statistical  properties  of  the  signals.  The 
details  of  the  simulation  approach  are  discussed  in  the  Appendix,  but 


highlights  of  the  approach  and  important  simulation  results  are  summarized  in 
this  section. 

The  measurement  approach  utilizes  a  laser  interferometric  system  for 
measurement  of  an  Interference  signal  which  is  sensitive  to  particle 
displacement.  Such  a  system  is  shown  schematically  in  Figure  1.  Although  the 
system  shown  is  physically  similar  to  the  familiar  Laser  Doppler  Velocimeter, 
important  differences  in  its  implementation  form  the  basis  of  the  present 
approach.  This  approach,  in  simple  terms,  is  to  measure  the  time-dependent 
interference  signal  arising  from  a  single  submicron  particle  (£0.1  pm 
diameter)  passing  through  the  measurement  volume  formed  by  the  intersection  of 
the  laser  beams  in  Figure  1  and  to  determine  the  particle  size  by  statistical 
analysis  of  the  signal.  To  understand  this  approach,  it  is  important  to 
appreciate  some  of  the  characteristic  aspects  of  the  Brownian  motion  of 
individual  particles. 

A  schematic  representation  of  the  time  dependent  motion  of  a  0.1  pm 
particle  is  shown  in  Figure  2  in  order  to  Indicate  the  time  scales 
characteristic  of  the  motion.  In  this  figure,  random  changes  in  particle 
velocity,  of  order  10"5  m«a-1,  occur  in  characteristic  times  of  10“14  seconds 
due  to  collisions  with  gas  molecules.  Given  sufficient  time  resolution  to 
measure  instantaneous  particle  velocities,  one  would  observe  a  Maxwellian 
velocity  distribution  if  a  sufficiently  large  sample  is  measured.  The 
characteristic  time  for  the  Maxwellian  distribution  to  evolve  is  the 
"relaxation  time"  t^^g*  (*  10”  7  s  in  Figure  2),  which  is  a  function  of 
particle  size  and  fluid  properties.  Clearly  if  one  wants  to  obtain  any 
information  on  the  Brownian  motion  of  an  individual  particle,  he  must  sample 
the  time  history  of  particle  motion  in  an  Interval  short  compared  to  the 
relaxation  time,  or  the  effects  will  be  averaged  out.  Conversely,  one  cannot 
monitor  the  Instantaneous  Brownian  velocity  of  the  particle,  because 
electronics  and  signal- to-noise  limitations  constrain  realizable  sampling 
Intervals  to  be  much  longer  than  the  10“ 14  s  characteristic  collision  time. 
The  measurement  strategy,  then,  is  to  sample  with  an  Interval  on  the  order  of 

s,  which  is  fast  compared  to  the  relaxation  time.  The  type  of  signal 


Figure  1.  Laser  Interferometric  System  for  Measurement 
of  Brownian  Motion. 
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Dependent  Motion  of  a  0.1  pm  Diameter  Particle 
r  at  300  K  and  Atmospheric  Pressure. 


measured  with  this  kind  of  time  resolution  will  be  a  superposition  of  many 
random  events,  but  the  time  dependence  of  this  signal  should  exhibit 
statistical  properties  representative  of  the  Brownian  motion  of  a  particle  of 
a  particular  size. 

A  second  Important  concept  is  mentioned  here  both  to  avert  possible 
misconceptions  about  the  ”LDV”  nature  of  the  measurement  and  to  underscore  the 
idea  that  the  signal  contains  time  (but  not  velocity)  information  which  must 
be  analyzed  in  a  statistical  sense.  This  concept  is  that  the  mean  excursion 
distance  for  the  particles  in  Brownian  motion  is  much  smaller  than  the 
wavelength  of  visible  light,  as  shown  in  Figure  3.  As  discussed  in  greater 
detail  in  the  Appendix,  the  consequence  of  this  characteristic  of  the 
measurement  is  that  the  magnitude  of  signal  variations  is  not  meaningful 
(because  the  motion  may  be  occurring  in  different  regions  of  the  LDV  fringe 
pattern).  The  frequency  characteristics  of  the  signal  are  significant, 
however,  in  that  they  can  be  related  to  the  statistics  of  particle  motion. 

The  identification  of  such  relationships  is  the  motivation  for  the  computer 
simulation  of  the  signal. 

The  computer  simulation  is  a  solution  to  the  equation  of  motion  for  a 
particle,  with  a  random  excitation  term  representing  random  collisions  with 
fluid  molecules.  A  Monte  Carlo  technique  is  used  to  simulate  the  velocity 
histories  of  such  particles  at  regular  sampling  Intervals  At.  Particle 
velocities  are  chosen  from  a  distribution  with  the  constraint  that  the 
distribution  is  Maxwellian  in  the  long  time  limit.  With  the  velocity  history 
simulated  in  this  fashion,  the  particle  trajectory  through  the  crossed  beam 
measurement  volume  of  the  Brownian  motion  sensor  can  be  determined  by 
integration  over  time  (a  mean  fluid  velocity  may  also  be  superimposed). 
Finally,  the  light  intensity  scattered  by  the  particle  from  each  of  the 
crossed  laser  beams  is  readily  calculated,  and  the  Interference  signal  at  the 
detector  may  be  obtained  as  a  function  of  time.  This  calculated  Interference 
signal  is  the  basic  result  of  the  simulation.  The  next  step  is  to  consider 
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methods  for  analysis  of  the  signal  In  order  to  relate  the  particle  size  to 
characteristics  of  the  signal. 

As  discussed  In  the  Appendix,  a  statistical  analysis  which  appears  to 
provide  a  suitable  correlation  Is  an  analysis  of  the  mean  time  between  local 
extrema  (first  derivative  zero  crossings)  In  the  signal.  As  shown  In 
Figure  4,  when  this  parameter  Is  plotted  against  a  parameter  depending  on  the 
sampling  Interval,  At,  and  the  particle  relaxation  time,  1/p,  a  straight  line 
correlation  Is  obtained.  (Curvature  away  from  the  linear  relationship  occurs 
when  the  sampling  Interval  becomes  comparable  to  the  relaxation  time).  Thus, 
from  a  measurement  of  a  statistical  parameter  (the  mean  time  between  extrema), 
the  particle  relaxation  time  and,  hence,  the  particle  mass  or  size  may  be 
determined.  This  approach  to  data  analysis  Is  the  approach  which  we  expect  to 
pursue  in  our  analyses  of  experimental  data. 

2.2  Experimental  Development 

The  experimental  development  during  the  first  year  of  the  program  has 
concentrated  on  the  design  and  assembly  of  a  laser  optical  system  suitable  for 
Interferometric  measurements  of  the  motion  of  small  particles  and  on  the 
design  of  a  data  acquisition  system  with  the  very  fast  data  rate  capabilities 
(“100  MHz)  required  by  the  measurement  technique. 

The  optical  system  which  has  been  assembled  for  this  study  Is  shown 
schematically  In  Figure  S.  A  Spec tra-Phy sics  Model  165  argon  Ion  laser  with  a 
single-line  power  of  approximately  2  watts  at  514.5  nm  is  used.  The  output  of 
this  laser  Is  expanded  and  collimated  to  a  beam  size  of  approximately  3  cm 
diameter  and  Is  then  split  with  a  50/50  cube  beamsplitter.  The  two  parallel 
beams  are  brought  to  a  common  focus  using  a  multi-element  lens  system.  The 
focused  laser  spot  size  (diagnostic  measurement  volume)  Is  approximately 
50  pm.  The  measurement  volume  is  Imaged  through  an  aperture  onto  a  fast 
photomultiplier  which  views  the  light  scattered  at  90°.  An  interference 
filter  is  used  to  reject  light  outside  a  1  nm  band  centered  at  514.5  nm. 
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Figure  4.  Correlation  Between  the  Mean  Time  Between  Signal  Extrema 
and  the  Particle  Relaxation  Time, P  . 


Figure  5.  Schematic  of  Experimental  Optical  System  for  the  Study 
of  Brownian  Motion  Measurement  Techniques. 


A  microcomputer-based  data  acquisition  system  with  very  fast  tine 
resolution  capability  has  been  assembled  for  this  project.  This  system  is  an 
IBM  Personal  Computer  XT  with  a  CAMAC  interface  (Transiac  Corp.,  Mountain 
View,  Calif.).  This  system  permits  considerable  flexibility  in  that  any 
CAMAC-compatlble  module  may  be  used.  Fast  time  resolution  in  data  acquisition 
is  achieved  with  the  appropriate  choice  of  CAMAC  hardware  in  combination  with 
fast  memories  for  temporary  data  storage.  The  contents  of  the  data  memories 
are  transferred  into  the  computer  memory  via  a  DMA  interface  at  data  rates  of 
approximately  400  kHz.  Large  volumes  of  data  may  then  be  stored  on  the 
system's  fixed  disk. 

Presently  we  have  set  up  and  tested  this  system  with  a  fast  amplifier  and 
transient  recorder,  enabling  acquisition  of  up  to  32,000  (8-bit)  data  points 
at  a  maximum  rate  of  100  MHz  (10  ns  sampling  interval).  We  anticipate  that 
this  recorder  will  be  used  in  our  Initial  experiments  for  testing  the  Brownian 
motion  measurement  concept.  Because  the  data  analysis  requires  statistical 
analysis  of  the  time  characteristics  of  the  signal,  we  may  ultimately  find 
that  a  time-domain  hardware  configuration  (e.g.,  for  time  distribution 
measurements  of  the  signal)  is  more  suitable.  The  CAMAC  system  provides  the 
flexibility  for  using  this  type  of  data  acquisition  hardware  without  redesign 
of  the  entire  data  acquisition  system. 

2.3  Research  Plan 

The  research  plan  for  the  second  year  of  the  program  emphasizes 
experimental  validation  of  the  measurement  concept.  Specifically,  we 
anticipate  performing  the  following  tasks: 

1.  We  will  initially  test  the  measurement  approach  by  measurement  of 
the  Brownian  motion  of  individual  submicron  particles  suspended  in 
liquids.  Commercially  available  particles  of  known  size  are  used  in 
these  measurements,  which  also  avoid  the  added  complications  of 
dispersing  the  particles  in  a  gas  stream.  These  tests,  in 
combination  with  the  Monte  Carlo  simulations,  should  establish 


whether  the  signal  fluctuations  exhibit  the  expected  behavior  for 
particles  of  different  size  and  should  provide  a  crucial  test  for 
the  measurement  theory. 

2.  If  the  tests  with  liquid  suspensions  conclude  with  favorable 

results,  we  will  proceed  with  measurements  of  submicron  particles 
suspended  In  gas  flows.  These  measurements  will  require  either 
dispersal  of  the  liquid-suspended  particles  Into  a  gas  stream  or 
perhaps  chemical  generation  of  small  particulates  from  the  gas 
phase.  Measurements  of  submicron  particles  suspended  In 
well-controlled  gas  flows  should  be  critical  In  establishing  the 
capabilities  and  limitations  of  the  technique  as  a  particle 
diagnostic  and  should  provide  the  foundation  required  for 
Investigations  of  diagnostic  applications  In  the  third  year  of  the 
program. 
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3.  PUBLICATIONS 


No  results  to  date  fron  this  program  have  been  published  in  technical 
journals.  We  anticipate  that  the  work  in  the  coming  year  on  experimental 
validation  of  the  Brownian  motion  measurement  approach  will  be  submitted  to  a 
journal,  e.g.,  Journal  of  Applied  Physics,  Applied  Optics,  Physics  of  Fluids, 
etc.,  for  publication. 
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The  Principal  Investigators  for  this  work  are  Dr.  Alan  Stanton,  who  is  a 
Principal  Research  Scientist  at  Aerodyne  Research,  Inc.,  and  Dr.  Ual  Cheng, 
Senior  Research  Scientist  (part-time)  at  Aerodyne  and  Assistant  Professor  of 
Mechanical  Engineering,  MIT.  Dr.  Stanton  has  concentrated  his  research 
efforts  to  date  for  this  program  on  the  development  of  optical  and  data 
acquisition  systems  for  the  experimental  program,  while  Dr.  Cheng  has 
concentrated  on  the  measurement  theory,  through  development  of  the  Monte  Carlo 
Brownian  motion  simulation  model.  In  addition  to  his  work  at  Aerodyne  on  this 
project,  Dr.  Cheng  has  supervised  the  work  of  undergraduate  and  graduate 
students  in  the  Department  of  Mechanical  Engineering,  MIT.  This  work  has  been 
in  support  of  the  Brownian  motion  simulation  model,  funded  by  a  subcontract  to 
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Three  presentations  were  made  at  AFOSR  Contractors'  Meetings  on  this  work 
during  the  first  year  of  the  contract.  These  presentations  were: 


A.C.  Stanton  and  W.K.  Cheng,  "Techniques  for  Submicron  Particle  Sizing", 
1983  AFOSR  Contractors'  Meeting  on  Air  Breathing  Combustion  Dynamics 
Research,  Scottsdale,  AZ,  September  19-22,  1983. 

W.K.  Cheng  and  A.C.  Stanton,  "Single  Particle  Sizing  by  Measurement  of 
Brownian  Motion",  1984  AFOSR  Research  Meeting  on  Diagnostics  of  Reacting 
Flows,  Yale  University,  March  21-22,  1984. 

A.C.  Stanton  and  W.K.  Cheng,  "Single  Particle  Sizing  by  Measurement  of 
Brownian  Motion",  1984  AFOSR/ONR  Contractors'  Meeting  on  Airbreathing 
Combustion  Research,  Pittsburgh,  PA,  June  20-21,  1984. 
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6.  NEW  DISCOVERIES  OR  INVENTIONS 


The  experimental  approach  under  study  in  this  program,  i.e., 
determination  of  particle  size  by  measurement  of  the  Brownian  motion  of 
individual  particles,  is  a  new  approach  which  has  not  previously  been 
demonstrated.  We  anticipate  that  this  approach,  if  successful,  will  permit 
the  development  of  research  instrumentation  for  the  study  of  the  growth  of 
submicron  particles,  based  on  new  concepts  in  measurement  and  data  analysis. 


APPENDIX 


1 .  Summary 

Approaches  to  signal  processing  for  the  Brownian  motion  sensor  are 
examined  by  simulating  the  signal  from  such  a  system  on  the  computer  using  a 
Monte  Carlo  method.  In  Section  2  the  nature  of  the  Brownian  motion  is 
discussed,  and  the  necessary  statistical  tools  to  generate  a  Monte  Carlo 
simulation  of  the  process  are  introduced.  Then  the  results  of  the  simulation 
and  the  qualitative  behavior  of  the  signal  from  a  laser  Doppler  system 
observing  the  particle  are  analyzed.  Because  the  amplitude  of  the  Brownian 
motion  is  smaller  than  the  fringe  spacing  of  the  probing  laser  beam,  the 
wavelength  of  the  laser  is  not  a  relevant  parameter,  and  only  temporal 
information  (and  hence,  not  velocity)  can  be  extracted.  In  Section  3,  the 
statistics  of  the  temporal  content  of  the  signal  are  examined.  We  find  that 
the  mean  period  of  the  signal  may  be  related  to  the  relaxation  time,  from 
which  the  mass,  or  size,  of  the  particle  may  be  calculated. 

2 .  Simulation  of  the  Brownian  Motion  Sensor  Signal 

A  numerical  simulation  of  the  signal  obtained  from  the  Brownian  motion 
sensor  is  used  to  study  the  nature  of  the  signal,  and  to  devise  an  algorithm 
for  extracting  the  information  pertinent  to  the  particle  size.  In  what 
follows,  the  method  of  simulating  the  signal  is  described  and  the  qualitative 
nature  of  the  signal  is  discussed.  Then  a  method  is  described  to  analyze  the 
signal. 

2.1  Nature  of  the  Brownian  Motion 

The  motion  of  a  particle  in  thermal  equilibrium  with  the  surrounding 
fluid  is  described  by  the  Langevln  equation,1 
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where  is  the  velocity  and  p  is  the  damping  coefficient.  The  random 
excitation  _A(  t)  has  the  property  that  the  solution  of  Eq.  (1),  which  is 
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has  the  Maxwellian  velocity  distribution  in  the  limit  of  large  t.  The  time 
scale  for  approaching  this  limit  is 


xrelax 


1/p  . 


(3) 


In  continuum  flow  (\mfp  «  d) ,  for  a  spherical  particle  of  diameter  d  and 
mass  m,  p  is  given  by 

p  -  3nn  d/m  ,  (A) 

where  ^  is  the  viscosity  of  the  fluid  medium.  For  the  free  flow  limit 
(\mfp  »  d),  the  value  of  p  is  given  by  the  Epstein  formula, 
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where  c  is  the  mean  thermal  speed  of  the  fluid  molecules  and  a  is  the 
accommodation  coefficient.  A  value  of  a  ■  0.5  is  used  in  the  simulation. 
Values  for  the  relaxation  time  are  shown  in  Figure  A-l. 

As  an  example,  we  consider  a  particle  of  0.2  ^m  diameter  with  a  specific 
gravity  of  2.  For  air  at  300  K  and  at  1  atmosphere  pressure,  the  relaxation 
time  is  Treiax  ■  150  ns.  To  obtain  information  about  the  velocity 


A-2 


RELAX 


distribution,  the  sampling  time  should  be  less  than  "trelax •  A  sampling  time 
At  of  10  ns  is  used  in  the  following  analysis.  Furthermore,  a  time  between 
samples  of  10  ns  is  assumed. 

The  time  scale  for  the  fluctuating  excitation  A(t)  is  of  the  order  of  the 
collision  time  xc  between  the  fluid  molecules  and  the  particle,  which  is 
~10~14  s.  During  each  collision,  the  particle  changes  its  velocity  by 

v8kTm  (6) 

m 

where  m  is  the  mass  of  a  fluid  molecule.  The  value  for  6V  is  ~10“5  m*s~^. 

Since  the  sampling  time  At  is  large  compared  to  the  collision  time  Tq, 
the  velocity  sampled  is  the  average  over  At  of  the  time  dependent  velocity 
given  by  Eq.  (2).  The  requirement  that  the  velocity  distribution  has  the 
asymptotic  limit  of  a  Maxwellian  distribution  requires  the  sampled  velocity 
Wq  at  nAt  to  have  a  probability  distribution  of,1 
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Therefore,  if  Wq  -  (un,vQ,wn),  and  if  the  Doppler  system  is  sensitive 
to  the  velocity  component  w,  the  probability  distribution  function  for  wn  is 
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2.2  Simulation  of  the  Brownian  Velocity 

Using  the  probability  distribution  of  Eq.  (8),  we  may  generate,  using  a 
pseudo-random  number  generator,  the  velocity  history  of  w  with  the  correct 
statistical  properties.  In  general.  If  R  Is  a  random  number  with  uniform 
distribution  In  the  Interval  [0,1  ],  and  p(x)  Is  the  probability  distribution 
of  a  variable  x,  R  may  be  mapped  into  x  by  equating  the  probability  of  finding 
x  m  [x1,x2]. 
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In  particular  let  Xj  «  -®  and  Rj  ■  0;  then 
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For  the  distribution  of  Eq.  (8),  It  Is  convenient  to  rewrite  the  expression  as 
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where 
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Then  Eq.  (10)  becomes 


—  /  exp  [-A2(w  -w  ,h)^1  dw  ■  R  ,  (14) 
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from  which  wn  may  be  solved  from  Che  random  number  R.  The  solution  of 
Eq.  (14)  Is 
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where 
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To  summarize,  the  algorithm  for  generating  the  Brownian  velocity  history 
works  as  follows.  The  initial  velocity  W}  Is  obtained  randomly  from  the 
Maxwellian  distribution 
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The  subsequent  velocities  wn,  with  sampling  time  At,  are  then  given  by 
Eqs.  (15)  and  (16). 
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2.3  The  Doppler  Signal 


Consider  a  crossed  laser  bean  system  with  beam  waist  r0  as  shown  in 
Plgure  1.  The  intensities  of  the  two  beams  scattered  by  the  particle  at  x 
(x,y,z)  (in  a  coordinate  system  centered  on  the  probe  volume)  are 
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Here  the  unit  propagation  vectors  are 
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Therefore  the  expressions  for  Ij  and  I2  become 
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The  mean  particle  trajectory  is 

*(t)  -  Xq  +  v  t  ,  (24) 

where  Xp  is  the  reference  point  and  v  is  the  mean  velocity.  Substitution  of 
x(t)  -  (x(t),  y(t),  z(t))  from  Eq.  (24)  into  Eqs.  (19)  and  (20)  would  give  the 
scattered  intensity  profile  as  the  particle  traverses  the  beam  volume. 

At  the  detector,  because  of  the  interference  between  the  two  beams,  the 
light  collected  is  a  modulated  Doppler  signal.  The  signal  may  be  written  as 

t 

S  ’  Ix  +  I2  +  cos  (Akz  •  /  v8dt  +  *)  ,  (25) 

where 

,  _  e 

Akz  -  (2  sin  f--)  ,  (26) 

and  the  z  velocity  component  contains  both  a  mean  and  random  velocity, 

vz  *  vz  +  w(t)  •  (27) 


The  phase  shift  4  is  a  constant  depending  on  the  reference  point  definition. 
The  power  received  at  the  detector  may  therefore  be  written  as 
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The  first  term  represents  the  absolute  intensity  function  which  Is 
proportional  to  the  beam  Intensity  Ic,  the  aperture  and  transmissivity 
function  H  of  the  receiving  optics,  and  the  scattering  cross-section 
integrated  over  the  appropriate  field  of  view.  The  functions  f1  and  f2  are 
defined  by  the  right  hand  sides  of  Eqs.  (19)  and  (20).  The  cosine  term 
represents  the  heterodyned  Doppler  signal. 

2.4  Simulation  Results  and  Discussion 

The  Brownian  velocity  of  a  0.1  pjn  diameter  spherical  particle  with  a 
specific  gravity  of  2  In  300  K  air  at  atmospheric  pressure  has  been  simulated 
using  the  Monte-Carlo  calculation  described  above.  The  velocity  history  Is 
shown  In  Figure  A-2  for  500  velocity  samples  with  a  sampling  time  of  10  ns 
each.  The  velocity  statistics  are  shown  in  the  histogram  In  Figure  A-3.  The 
result  agrees  with  the  Maxwellian  distribution  which  Is  also  plotted  in  the 
figure. 

Calculations  were  then  carried  out  to  examine  the  behavior  of  the  Doppler 
signal,  which  Is  described  by  Eq.  (28).  The  term  of  Interest  Is 

S  “  ltfl  +  f2  +  2^1*2  COS  (TT  *  2  8ln  T’  tV  +  f  +  ♦>}  *  (29) 

—  CD 

with  f}  and  f2  defined  by  the  right  hand  sides  of  Eqs.  (19)  and  (20).  To 
obtain  the  most  favorable  signal,  the  particle  Is  assumed  to  have  mean 
velocity  In  the  y  direction  only,  and  is  traversing  the  center  of  the  sample 
volume  along  the  y  axis.  The  reference  phase  value  is  set  arbitrarily  to 
(-0.5x)  which  corresponds  to  aligning  the  y  axis  along  the  "boundary"  of  the 
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PARTICLE  DIAMETER 


Figure  A-2.  Simulated  Brownian  Velocity  History  of  a  0.1  ym 

Diameter  Particle  (Specific  Gravity  =  2)  in  Air  at 


Comparison  of  the  Velocity  Distribution  Obtained  from 
the  Monte  Carlo  Simulation  with  a  Maxwellian  Distribution 


bright  and  dark  fringes  In  the  beam  sample  volume.  The  laser  beam  waist  rQ 
is  50  pm  and  the  beam  crossing  angle  is  set  at  64  -  10s.  These  values  are 
typical  of  a  practical  system. 

Under  the  above  conditions,  the  simulated  Doppler  signal  for  a  0.1  pm 
particle  with  vy  «  1  m*s-1  traversing  the  center  part  of  the  beam  volume, 
starting  from  x  -  (0.,  -10  pm,  0.),  is  shown  in  Figure  A-4(a).  The  horizontal 
axis  is  the  number  of  time  steps,  with  each  time  step  equal  to  10  ns.  The 
signal  fluctuation  due  to  the  Brownian  motion  is  clearly  detectable. 

The  simulation  is  repeated,  under  the  same  conditions,  with  a  0.01  pm 
diameter  particle.  The  simulated  signal  for  the  smaller  particle  is  shown  in 
Figure  A-4(b) .  The  fluctuation  is  of  a  much  higher  level  than  those 
exhibited  by  the  larger  diameter  particle.  Furthermore,  because  of  the 
behavior  of  the  cosine  function,  the  modulation  is  reduced  when  the  argument 
of  the  cosine  is  close  to  ±  an,  a  -  0,1,2,....  For  the  signal  in  Figure 
A-4(b),  the  particle  is  very  close  to  the  center  of  the  beam  volume  as  the 
absolute  intensity  is  approximately  constant.  Therefore,  the  signal  has  a 
magnitude  of 


S  ~  1  +  cos  {^2  .  2  sin  [v^t  +  /  w(t)dt]  +  .  (30) 

Thus,  at  regions  where  S  ~1  or  2,  the  magnitude  of  the  modulation  is 
reduced.  This  effect  is  evident  in  Figure  A-4(b). 

The  signal  modulation  in  Eq.  (29)  is  a  result  of  the  cosine  term,  which 
is  rewritten  here  as 

modulation  ~  cos(k'z  +  $)  (31) 
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Figure  A-4.  Simulated  Time  Dependent  Interference  Signals 
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(b)  0.01  pm  Diameter  Particle 
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The  root-mean-square  value  <z>rm8  of  the  Brownian  displacement  Is  shown  In 
Figure  A-5.  For  particles  of  100  to  1000  A  diameter,  <zm>/X  “  0.1. 

Therefore  the  rms  displacement  Is  much  less  than  the  fringe  spacing  k,_1.  As 
a  result,  the  length  scale  X  does  not  play  any  role  In  the  determination  of 
the  Brownian  velocity  w.  The  signal,  therefore,  only  contains  Information 
about  the  time  scale  and  not  any  length  scale,  and  the  determination  of  a 
meaningful  velocity  Is,  In  principle,  not  possible. 

The  essence  of  the  above  discussion  Is  illustrated  In  Figure  A-6,  which 
exhibits  the  transfer  function  of  the  detection  system  [the  right  hand  size  of 
Eq.  (29)].  For  typical  Brownian  motion  modulation  (AB),  the  signal  modulation 
is  (A'B').  The  signal  modulation  is  drastically  reduced  where  the  slope  of 
the  transfer  function  vanishes  (CD).  The  amplitude  of  the  modulation  is  a 
function  of  the  position  of  the  particle  relative  to  the  fringe  pattern. 

Since  this  Is  not  known  a  priori,  the  amplitude  modulation,  therefore,  does 
not  provide  useful  information  on  the  nature  of  the  Brownian  motion. 

When  there  is  a  fluid  velocity  component  v8  perpendicular  to  the  fringe 
pattern,  the  particle  position  will  be  represented  by  the  trajectory  (EF), 


lion  Distance,  Compared  with  the  Wavelength  of  Visible 
Particles  (Specific  Gravity  =  2)  in  Brownian  Motion  in 
=  300  K,  p  =  1  Atmosphere. 


Modulation  of  the  Doppler  Signal  as  a  Function  of 
Particle  Position. 


with  Che  Brownian  modulation  such  as  (AB)  superimposed  on  Cop.  The  signal 
will  exhibit  full  modulation  (E'F‘)  (regular  Doppler  signal)  with  the  Brownian 
motion  modulation  (A'B1)  superimposed  on  top.  When  vz  Is  large  compared  to 
the  Brownian  velocity  w.  It  may  be  difficult  to  recover  the  Brownian  motion 
modulation  from  the  regular  Doppler  signal.  This  effect  Is  Illustrated  In 
Figures  A-7(a)  and  A-7(b)  In  which  velocity  7Z  of  0.1  and  0.3  m«s-1 
respectively  are  superimposed  onto  the  Brownian  velocity.  In  the  latter  case, 
the  Brownian  motion  signal  is  almost  completely  masked  and  the  signal 
resembles  the  usual  Doppler  velocimeter  signal.  For  small  particles  with  a 
higher  Brownian  velocity,  good  quality  signal  may  still  be  obtained  with  a 
significant  vz.  This  is  illustrated  in  Figure  A-8  by  the  simulation  of  the 
signal  from  a  0.01  nm  particle  with  "vz  -  0.3  m^s"1.  The  mean  Brownian 
velocity  w,  where 


Is  approximately  0.1  m^s”1  for  a  0.1  ^m  diameter  particle  at  a  gas  temperature 
of  300  K.  Thus  one  expects  that  for  particles  In  the  0.01  to  0.1  ^m  diameter 
range,  reasonable  signal  quality  may  be  obtained  If  "vz  Is  less  than 
0.1  ra»8_1. 

The  above  discussion  indicates  that  the  strategy  for  analyzing 
characteristics  of  the  Brownian  motion  lies  in  analysis  of  the  time 
Information.  The  fundamental  time  scale  of  the  fluctuation  is  the  relaxation 
time  1/0.  The  information  on  the  particle  size  may  be  calculated  from  6  using 
Eqs.  (4)  and  (5).  Furthermore,  the  signal  must  be  filtered  to  remove  the 
regular  Doppler  component  (vz)  for  proper  Interpretation. 


3.  Processing  of  the  Brownian  Motion  Sensor  Signal 
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Figure  A-7.  Simulated  Brownian  Motion  Signals  for  a  0.1  ym  Diameter  Particle  with 
Fluid  Velocity  Component  ^  Superimposed 
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particle  due  to  Its  mean  velocity),  the  latter  may  be  filtered  out  easily.  In 
the  following  discussion  we  shall  assume  that  the  signal  has  already  been 
filtered  to  contain  the  Brownian  Information  only.  The  appropriate  statistics 
of  the  temporal  information  signal  are  discussed. 

The  time  scale  of  the  fluctuating  signal  may  be  analyzed  in  terms  of  the 
time  between  zero  crossings  (assuming  that  the  mean  has  already  been  removed 
by  subtraction  of  a  moving  average),  or  in  terms  of  the  time  between  extrema 
of  the  signal  (zero  crossing  of  the  first  derivative).  The  particular 
implementation  depends  on  the  hardware  configuration.  Nevertheless  the 
principle  remains  the  same.  Ve  shall  use  the  time  between  extrema  method  in 
the  following  analysis.  Also  the  analysis  considers  only  the  mean  time 
between  extrema,  and  does  not  attempt  analyze  the  full  spectral  content  of  the 
signal,  so  that  possible  aliasing  is  not  an  Issue.  (The  method  amounts  to 
doing  statistical  sampling  rather  than  fixed  Interval  sampling,  so  on  the 
average,  the  higher  frequency  content  la  "represented.") 

The  signal  processing  Is  simulated  on  the  computer  using  the  Monte  Carlo 
technique  described  previously.  The  Brownian  motion  of  a  particle  with 
specific  gravity  of  two,  embedded  In  room  temperature  air,  is  analyzed.  The 
mean  time  between  extrema  t  as  a  function  of  the  particle  diameter  d  Is  shown 
In  Figure  A-9.  The  values  of  T  do  not  seem  to  follow  a  simple  linear 
relationship  with  d,  notwithstanding  the  fact  that  the  only  time  scale 
Involved  is  the  relaxation  time  Treiax  ■  1/p.  With  p  proportional  to  1/d  In 
the  free  flow  regime,  the  relaxation  time  should  be  proportional  to  d. 
Furthermore,  depending  on  the  sampling  interval  At,  the  mean  time  "t  Is 
different.  The  relevant  parameter  should  therefore  clearly  involve  At.  This 
Issue  Is  resolved  when  the  data  are  replotted  against  -^At/p  (with  the 
dimension  of  time)  (Figure  A-10).  The  data  points  collapse  on  to  a  single 
straight  line  passing  through  the  origin.  The  relationship  between  sampling 
and  adequate  statistics  is  also  illustrated,  as  at  large  values  of  ^At/p,  a 
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Figure  A-9.  Mean  Time  Between  Extrema  in  Brownian  Motion 

Interference  Signal,  for  Two  Sampling  Intervals. 


b 


A-21 


1  .  ' 


r*  ** 


Figure  A-10.  Correlation  Between  the  Mean  Time 
and  the  Particle  Relaxation  Time, 


long  string  of  data  Is  required  to  produce  a  good  estimate  of  t.  For  small 
values  of  "^At/p,  the  sampling  Interval  of  At  ■  10  ns  is  not  adequate  to 
resolve  the  Brownian  motion  time  scale,  and  as  a  result,  t  Is  overestimated 
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